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A state in quantum mechanics is defined as a positive operator of norm 1. For finitc 
systems, this may be thought of as a positive matrix of trace 1. This constraint of 
positivity imposes severe restrictions on the allowed evolution of such a state. From 
the mathematical viewpoint, we describe the two forms of Standard dynamical 
equations - global (Kraus) and local (Lindblad) - and show how each of these 
gives rise to a semi-group description of the evolution. We then look at specific 
examples from atòmic systems, involving 3-level systems for simplicity, and show 
how these mathematical constraints give rise to non-intuitive physical phcnomena, 
rcminiscent of Bohm-Aharonov effects. In particular, we show that for a multi- 
level atòmic system it is generally impossible to isolate the levels, and this leads 
to observable effects on the population relaxation and decoherence. 



1. Introduction 

The Standard description of a quantum state suitable for an open system 
is by means of a density matrix p, a positive matrix of trace 1. For a 
hamiltonian (non-dissipative) system one obtains a unitary evolution of 
the state. For a non-dissipative system the time evolution of the density 
matrix p(t) with p(to) = po is governed by 

p(t) = U(t) Po U(t)\ (1) 
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where U(t) is the time-evolution operator satisfying the Schrodinger equa- 
tion 

ihjU{t) = HU(t), U(0) = I, (2) 

where / is the identity operator. The state p(t) equivalcntly satisfies the 
quantum Liouville equation 

ihj t p{t) = [H, p(t)\ = Hp(t) - p(t)H. (3) 

H is the total Hamiltonian of the system. (In the context of Quantum 
Control theory, we may assume that H = H(f) depends on a set of control 
fields f m : 

M 

H(f)=H +J2fm(t)H m , (4) 

m— 1 

where Hq is the internal Hamiltonian and H m is the interaction Hamiltonian 
for the field f m for 1 < m < M.) The advantage of the Liouville equation 
(3) over the unitary evolution equation (1) is that it can easily be adapted 
for dissipative systems by adding a dissipation (super-)operator Lr>[p{t)]: 

ihp(t) = [H,p(t)]+ihL D [p(t)}. (5) 

In general, uncontrollable interactions of the system with its environ- 
ment lead to two types of dissipation: phasc decoherence (dephasing) and 
population relaxation (decay). The former occurs when the interaction with 
the enviroment destroys the phase correlations between states, which leads 
to a decay of the off-diagonal elements of the density matrix: 

1 

Pkn{t) = --([H,p(t)]) kn - T kn p kn {t) (6) 

where Tfc„ (for k ^ n) is the dephasing rate between |A:) and \n). The latter 
happcns, for instance, when a quantum particle in state \n) spontancously 
cmits a photon and decays to another quantum state |fc), which changes 
the populations according to 

Pnn(t) = --{[H,p{t)]) nn + ^ hnkPkk(t) - JknPnn(t)] (7) 

where jknPnn is the population loss for level |n) due to transitions |n) — * 
|fc), and jnkPkk is the population gain caused by transitions |fc) — > \n). 
The population relaxation rate 7fc„ is dctcrmined by the lifetime of the 
state \n), and for múltiple decay pathways, the relative probability for the 
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transition \n) — > Phase decoherence and population rclaxation lead 
to a dissipation superoperator (represented by an N 2 x N 2 matrix) whose 
non-zero elements are 

{L D )kn,kn = -^kn k ^ n 

(L D )nn,kk = +lnk k^íl (8) 

{L[^) nnnn ^In^k ^kn 

where Ykn and 7fe„ are positive numbers, with T^ n symmetric in its indices. 
The N 2 x N 2 matrix superoperator L d may be thought of as acting on the 
iV 2 -vector V obtained from p by 



V[(i-l)N+j] = Pij- 

The resulting vector equation is 

V = LV = (L H + L D )V 



(9) 



(10) 



where Lh is the anti-hermitian matrix derived from the hamiltonian H. 

The vàlues of the rclaxation and dcphasing parameters may be deter- 
mined by experiment, or simply chosen to supply a model for the dissipa- 
tion phenomenon. But they may not be chosen arbitrarily; the condition of 
positivity for the state p imposes constraints on their vàlues, as does their 
dcduction from rigorous theory. We illustrate this by demonstrating the 
constraint for a two-level system. 



2. Two-level systems 
2.1. Unitary evolution 

The general hamiltonian for a two-level system is given, up to an additive 
constant, by 



"1 


0" 


+ fx 


'0 1" 


+ fy 


'0 -i' 





-1 


.i o_ 


i 



H = w 

which we assume here to be time-independent. 

This corresponds to the (superoperator form) Lh, where 

"0 Í(fx+Ífy) Í{-fx+Ífy) 

-i (~fx + ify) -2ÍW i (~f x + if y ) 

-i (fx + ify) 2 ÍW Í(fx+ ify) 

.0 -Üfx+ify) -Í(~fx+ify) 



(11) 



L H = 



(12) 



February 1, 2008 10:59 WSPC/Trim Size: 9in x 6in for Procccdings 



DGMTP23v2 



4 



Note the useful rule for obtaining the equivalent iV 2 x N 2 matrix action 

ApB A ® B T V. 

The corresponding evolution equation for the 4-vector V corresponding 
to the state p is 



V = L H V. 



(13) 



This is equivalent to Eq.(l), which clearly preserves the trace of p, and also 
its positivity, using the dcfinition of a positive matrix as one of the form 
MM^ . (Of course this result is true in general.) 



2.2. Pure dissipation 

The dissipation (super-)operator is 



"-721 








712 





-r 














-r 





. 721 








-712 



The corresponding evolution equation 

V = L D V. 

has solution 

V(t) = exp(L D t)V(0) 
which corresponds to a value of the state p(t) 



Pll(7l2+72l-E)+7l2P22(l--E) 

721+712 

e~ tr p 21 



72lPll(l--E)+P22(72l+7l2-E) 

721+712 



(14) 



(15) 
(16) 

(17) 



where E = e t (T 21 +T 12 ) j for which it may readily be checked that Trp(í) = 
P11 + P22 = 1- Additionally, detp(í) is given by 

v 2 



PiiP22e- t(721+7l2) - (e- 2tr 



) P12P21 + 2 



P11712P22721 (1 - e-^+^y 



(721 +712)^ 



which is clearly positive for all t when 

2r > 712 + 721 



(18) 
(19) 
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since detp(i) > e _í ' 721+712 ' det p(0) > 0. Conversely, when the condition 
Eq.(19) is violated, it is easy to display examples for which the evolution 
does not produce a state. For example, for a pure state, which satisfies 
P11P22 - P12P21 = 0, choosing 712 > 2r, 721 = 0, Eq.(18) is clearly negativc. 

2.3. General dissipation 

When the hamiltonian matrix Lh and the dissipation matrix Lr> com- 
mute, the conclusions of the previous two subsections produce the same 
constraint for the solution of Eq.(10). In the general case these matrices do 
not commute; they do however generate a local semi-dircct group. More 
accurately the Lie àlgebra is locally a semi-direct sum a , which then gener- 
ates a semi-group. In this case also, general theory, which we discuss in the 
next section, shows that the trace and determinant conditions of Eq.(19) 
remain unchanged. 

3. Rigorous formulations 

3.1. Kraus formalism and semi-groups 

The global form of the evolution equation Eq.(l) in the presence of dissi- 
pation is due to Kraus 2 . The evolution of the state p is given by 



Equation(20) and the condition Eq.(21) clearly guarantee both positivity 
and unit trace. 

Further, though less obviously, this system implies the existence of a 
semi-group description of the evolution. For if we consider the set G whose 
elements are the sets {wí} satisfying Eq.(21), then if g = {vüí} and g' = {w-} 
are two elements of G, then so too is gg', where the product is taken in the 
sense of set multiplication. Although closed under composition, the only 
elements of G which possess inverses are the singleton sets {U}, where U 
is unitary. 



a In the present two-level case, the local Lie àlgebra is the 12-element gl(3, R) © R , and 
in general gl(N 2 — 1, R) © R N _1 , as discussed in 1 . 




(20) 



with 




(21) 
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3.2. Lindblad formalism 

In so far as the Kraus formalism provides an analogue of the unitary evo- 
lution equation Eq.(l), the Lindblad 3 formalism gives an analogue of the 
Schroedinger equation Eq.(3): 

p(t) = L[p{t)]p{t) 

= -i[H,p{t)] + ±£ (\y k p(t),V*] + [V klP {t)vl\) (22) 

k 

where the \% arc N x N matrices, but othcrwise arbitrary b . It may be 
proved that the dissipation superoperator arising from Eq.(22) has neg- 
ative eigenvalues. Since the evolution dynamics arises from exponentiation 
of Ljjt it follows that operators expfX/jí) in the theory will become un- 
bounded for arbitrary negative t. This means that not all operators will 
have inverses and implies a semi-group character to the evolution, as in the 
Kraus formalism. 

3.3. 2x2 Lindblad example 

Choosing four independent complex F-matrices 





"oi 0" 




"0 a 2 ~ 




"0 0" 




"0 0" 


Vi = 


0_ 


v 2 = 


0_ 


v 3 = 


_a 3 0_ 


v 4 = 


a 4 



we obtain for the dissipation superoperator Lu 



"-l«3l 2 








|a 2 | 





-1/2,4 














-1/2A 





. M 2 








-|«2 



where A = \a\\ 2 + | «2 1 2 + | a. 3 1 2 + | a.4 1 2 , which on comparison with Eq.(14) 
gives, defining f = ^(|ai| 2 + |a 4 | 2 ) 

72i = |a3| 2 , 7i2 = |«2| 2 , r = f+ ^ (712 +721) 

whence the constraint Eq.(19). Note that (712 + 72i)/2 is the phase deco- 
herence forced by population relaxation and T is the contribution of pure 
dcphasing. 



'We may also choose an arbitrary number of matrices Vfc . 
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3.4. General N X N Lindblad case 

A convenient choice for the 14 matrices may be made by defining 

V [i,j] = a [i,j] E H 

where Eij is the Standard basis for N x N matrices, with (_Ejj) a/ g = 5i a 5jp 
and we use the index notation = (i — 1)N + j. The relaxation and 
decoherence parameters are defined by 



Jij = 


i i 2 c 


Ï3) 


f — 


1 ,2 

2 M a [j,í] | "+ 


\ a U,j]\ 2 ) (i^i) 


r — 


lj2(\ a [k,i 
z i. 


+ \ a [k,j] ) (i^j) 



4. Bohm-Aharonov type effects 

What we mean by Bohm-Aharonov type effects in the title of this note, and 
of this section, is the impossibility of isolation of quantum subsystems. We 
illustrate this type of effect by considering the use of a two-level atòmic 
system as, say, a qubit, when this is a subsystem of a multi-level system. 

We consider the case of pure dissipation as discussed in subscction 2.2. 
Choosing vàlues 721 = 0, 712 =7, T = ^7, which satisfy the constraint 
Eq.(19), the state evolution is given by 



P(t) 



pu + P22 (1 



-*7 



,-1/2 t 7 



P21 



where the initial state is 



p(0) = 



Pu P12 
P21 P22 



„-l/2í-y - 



P22& 



-í 7 



(24) 



We now assume that our two-level system is embedded in a three level 
system, so that the state's evolution is given by 



(25) 



Now consider three examples for the state evolution. In all cases we start off 
with a pure state, in the first case with the third level not being populated. 





' P11 + P22 (1 - e t7 ) 


e-V2*7 pi2 


Pl3 


p(t) = 


e- 1/2t >2i 


P22e _t7 


P23 




_P31 


P32 


P33. 
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4.1. Unpopulated third level 

Assume an initial pure state represented by the 3- vector v = [1/ y/2, 1 / y/2, 0] 
corresponding to the density matrix 



1 1 
1 1 




Assuming that the third level is unaffected, the state evolution is given by 
(measuring t in units of I/7) 

'l-l/2e-' l/2e-V 2t 0" 

172è- 1 / 24 l/2e- 4 



P(t) = 















(26) 



In this case the naive picture of the evolution is justified, as the third level 
remains unpopulated, the eigenvalues remain positive (> 0), and the extra 
levels are not affected by the dissipative dynamics. The third level plays no 
role in the evolution. However, in general an upper level will not be totally 
unpopulated; and in this case the constraints play a role. 



4.2. Equally populated third level 

We take the initial pure state vector to be 

v = [Í/VS, l/y/3, 1/V3] 

giving the evolution 

- 2 _ e -t e -V2í 

P(t) 



-1/2Í 



(27) 



(28) 



As in subsection 4.1 we have assumed that the third levels are not affected 
by the dissipative dynamics. However, a numerical calculation shows that 
the eigenvalues of p(t) are not all positive; therefore the assumed evolution 
does not give a state, and so the naive assumption that the other levels 
remain unaffected is false. 



4.3. Pure dephasing 

Population relaxation is not the only source of constraints on the decoher- 
ence rates for N > 2. Even if there is no population relaxation at all, i.c., 
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7fe„ = for all k,n, and the system experiences only pure dephasing, we 
cannot choose the decoherence rates Tkn arbitrarily. For example, sctting 
Ti2 7^ and T 2 3 = Ti 3 = for our three-level system gives 



p(t) = 



Pu 
-r i2 í 

P31 



p21 



Pl2 Pl3 
P22 P23 
P32 P33 . 



(29) 



Choosing p(0) as in Eq. (27) we again obtain a density operator p{t) with 
negative eigenvalues, as a simple calculation will reveal. This shows that 
there must be additional constraints on the decoherence rates to ensure 
that the state of the system remains physical. 



5. Conclusions 

We have shown that it is impossible to isolate a two-level system from a 
multi-level system in the sense of assuming that the other levels will not be 
affected by relaxation and decoherence effects in the "isolated" system. A 
more general treatment of the effects noted here may be found elsewhere 4 ; 
in that paper the constraints are explicitly described for some multilevel 
systems, and the effects of these constraints are discussed. 
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